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ABSTRACT 

Aims. We present a sample of candidate quasars selected using the A'X-technique. The data cover 0.68 deg^ of the X -ray Multi-Mirror 
(XMM) Large-Scale Structure (LSS) survey area where overlapping multi-wavelength imaging data permits an investigation of the 
physical nature of selected sources. 

Methods. The KX method identifies quasars on the basis of their optical (R and z') to near-infrared (A',) photometry and point-like 
morphology. We combine these data with optical {u*,g'r', i', z') and mid-infrared (3.6 - 24 jim) wavebands to reconstruct the spectral 
energy distributions (SEDs) of candidate quasars. 

Results. Of 93 sources selected as candidate quasars by the KX method, 25 are classified as quasars by the subsequent SED analysis. 
Spectroscopic observations are available for 12/25 of these sources and confirm the quasar hypothesis in each case. Even more, 90% 
of the SED-classified quasars show X-ray emission, a property not shared by any of the false candidates in the KX-selected sample. 
Applying a photometric redshift analysis to the sources without spectroscopy indicates that the 25 sources classified as quasars occupy 
the interval 0.7 < z < 2.5. The remaining 68/93 sources are classified as stars and unresolved galaxies. 

Key words. Photometry - Quasars ; general - Surveys 



1. Introduction 

Despite the enormous progress that wide-field imaging surveys 
have provided in understanding quasar populations over the past 
decade, concerns remain over the extent to which observational 
biases contribute to a partial view of the quasar phenomenon. 
Consequently, multi-wavelength observations are important not 
only to reproduce the spectral energy distribution (SED) of the 
objects under study, but also to provide as complete a census as is 
practical of the quasar population. Nevertheless, selecting source 
populations using different energy bands encompasses a risk of 
confusion between different types of objects or different physical 
processes. Thus, studying the selection effects is an essential step 
to understanding the properties of the parent population. 

A classic method to identify type-1 quasars in optical imag- 
ing surveys is known as the f/V-excess (UVX; Sandage 1965 
Schmidt & Green||1983| . The UVX technique exploits the fact 
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that quasars display an emission excess in blue wavebands com- 
pared to main sequence, late type Galactic stars, and therefore 
occupy a distinct (i.e. bluer) locus in a color-color diagram with 
respect to stars. 

Based on a study of 323 radio-loud sources selected from the 
Parkes catalog ( Wright & Otrupcek|1990| l, approximately 70% 
of which were identified spectroscopically as quasars, 'Webster] 
|et al. ( 1995 1 claimed that up to 80% of quasars could be missed 
in optical surveys. Their suggestion was based on the Bj - K 
color index of these radio-selected sources, which showed a 
much larger scatter (1 < Bj-K < compared to the typical val- 
ues for optically selected quasars. Various interpretations were 
suggested, such as (a) intrinsic light absorption, (b) absorption 
from dust along the line of sight or (c) differences in the internal 
processes in the active nucleus. Cases (a) and (b) would result 
in "reddened" quasars, whilst objects in category (c) would be 
considered as intrinsically "red" (in contrast to the "classical", 
blue type-1 quasars). 

However, a number of studies have questioned the dusty 
quasar hypothesis. |Masci et al.| ( [T998| rejected the dust hypoth- 
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esis and |Whiting et al.| ( |2001| l suggesting that the B j - K spread 
could be accounted for by synchrotron emission. 'Benn et al. 



( |1998| l claimed that the Bj - K scatter could be attributed to 
observational effects such as variability, underestimated photo- 
metric uncertainties, etc. Althoug h jFrancis et al. j2000J support 
the conclusions of [Webster et al.| ( |1995| ), [Richards et al.| ( |2003"l l 
demonstrate that, for the majority of the SDSS quasars, redder 
colors could be explained by intrinsic processes in the AGN and 
that the number of missed red and reddened quasars in surveys 
such as the SDSS (different from typical optical surveys because 
of the filter selection) comes down to "only" 15%. Brown et al. 
(|2006 ) found similar results, suggesting that red type-1 quasars 
correspond to about 20% of the type-1 population. 

With the advent of large field, near infrared (NIR) surveys. 
Warren et al.| 1*2000 ) suggested a new colour excess criterion, de- 
signed to be less prone to the problem of dust reddening and ex- 
tinction. Their method is the NIR analog of UVX: as the quasar 
SED follows a power-law, it displays a flux excess at NIR wave- 
lengths compared to the Rayleigh-Jeans tail observed in stars. As 
a result quasars should also occupy a distinct region of the NIR 
color-color diagram, as they do in the UBV color-color plane. 
The so-called A'-excess criterion, hereafter KX, would serve 
as an alternative to the UVX for selecting quasars at redshifts 
z > 2.2. At these redshifts, the UV/optical colors of quasars be- 
come virtually indistinguishable from those of late type stars. 

Despite being a prominsing technique, relatively few KX- 
selected quasar samples have appeared in the literature - mainly 
due to the considerable effort required to complete a large 
area NIR imaging survey. The conclusions drawn from the 
/TX-selected quasar populations published to date are rather 
limited as various factors such as the small area of the sur- 
vey (|Croom et al. 2001), the relative shallow magnitudes sam- 
pled (Barkhouse & Hall 2001) or the inhomogeneity and in- 
completeness of the survey (Sharp et al. 2002| l have hampered 
an assessment of the potential of the method. The first detailed 
study of the KX method was given by Jurek et al.' (2007 ) while 
more recently [S mall et al. (2008 ) also published a work on 
/TX-selected quasars, raising the problem of contamination by 
foreground compact galaxies and their impact on spectroscopic 
follow-up observations. 

In this paper we present an analysis of a sample of 93 KX- 
selected quasar candidates selected from 0.68 deg.^ of Rz'K, 
imaging data located in the X-ray Multi-Mirror (XMM) Large 
Scale Structure (LSS) survey field ( [Pierre et al.| 2004). The 
XMM LSS survey is an X-ray imaging survey covering ap- 
proximately 10 deg^ to an approximate point-source flux limit of 
8 X lO^'^ergs^'cm"^ in the [0.5-2] keV energy band. The XMM 
LSS is associated with a number of imaging data sets at dif- 
ferent wavebands: the XMM LSS contains the XMM Medium 
Deep Survey (XMDS), a deeper X-ray imaging component 
over 2 deg.^ in addition to the optical Canada-France-Hawaii 
Telescope Legacy Survey (CFHTLS) and the Spitzer Wide area 
IR Extragalactic (SWIRE) survey. These additional data provide 
a panchromatic view of the quasar candidates and permit a de- 
tailed discussion of their physical nature. 

The paper is organized as follows: Section 2 presents the 
optical and NIR data used in the KX technique. Section 3 de- 
scribes the multi-wavelength data available in the field. Section 
4 describes the KX selection criteria and discusses the physi- 
cal nature of the selected sources using their multi-wavelength 
properties. In Section 5 we draw some conclusions on the color 
properties of the confirmed quasars selected using KX and com- 
ment on the overall efficiency of the technique. We also consider 



some of the challenges to be overcome in applying the KX tech- 
nique to large data sets. 



2. Optical and near-infrared observations and data 
reduction 

Two data sets were combined to produce a single optical 
plus near-infrared catalog upon which to conduct a search for 
quasars using the KX technique. Optical R,z' data were ob- 
tained at the CeiTO Tololo Interamerican Observatory (CTIO) 
and near-infrared /Ts-band data were obtained at Las Campanas 
Observatory (LCO). The two data sets were reduced indepen- 
dently and subsequently matched to form a single catalog. In 
this section we briefly describe the observations, data reduction 
and general properties of these photometric data sets. 



2.1. Near-infrared observations 

Near-infrared observations were performed during the period 
25 - 28 October 2002 using the 2.5 m (100 inch) Du Pont tele- 
scope, at Las Campanas Observatory. The XMDS field was ob- 
served in the A'j-band, using the Wide Field Infrared Camera 
(WIRC, [Persson et aL][2002] l. WIRC is a NIR (1.0 - 2.5 yum) 
camera consisting of four 1024^ Rockwell HAWAII arrays. Each 
array covers an area 3.4 x 3.4 arcmin^ with a scale of 0.2"/pixel. 
The four detectors are equally spaced with respect to the center 
of the configuration, with a gap between them of about 3.05'. 

Each telescope pointing consisted of a regular grid of nine 
dithering pointings with a 20" offset. A total of 4 x 28s expo- 
sures were obtained at each dithered location and each combined 
image displays a total integration time of 9 x 4 x 28s = 1008s 
in the central area. The nights of Oct. 25 and Oct. 27 were of 
mediocre photometric quality, while the nights of Oct. 26 and 
Oct. 28 were photometric. The typical seeing conditions during 
the run were between 1.1 - 1.5" FWHM. Several hundred im- 
ages were obtained and then combined to generate a mosaic of 
0.8 degl 



2.2. Optical observations 

The full XMDS region was imaged in the R- and z'-bands using 
the Mosaic 11 CCD imager mounted on the 4 m Blanco telescope 
at CTIO observatory. The Mosaic II camera consists of eight 
2048 X 4096 CCD arrays. The scale at the center of the camera 
is 0.27"/pixel, covering in total 36' x 36' on the sky. The ob- 
servations were obtained in November 2001 during photometric 
conditions with an average seeing of FWHM 1.3". Standard 
procedures were followed for the reduction of the obtained im- 
ages. 

The CTIO individual frames were astrometrically calibrated 
using the USNO-2 catalog (Monet et al. 2003) in order to gen- 
erate mosaic images of the observed area. The final astromet- 
ric precision reached was on the order of 0.3". Source extrac- 
tion and photometry were performed using SExtractor ^Bertin[ 
|& Arnouts 1996) with independent catalogs generated for each 
band. Detailed information regarding the data processing of the 
CTIO data can be found in [Andreon etaL] ( |20^ . 
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2.3. Near-infrared data processing 

Near-infrared data were proce ssed using a set o f scripts based 
upon the IRAPQand PHIIRS (" Hall et al.||T998| packages. The 
reduction strategy is based upon previous near-infrared surveys 
reported by Chen et al. (20^ and |Labbe et al. (2003 ). A com- 
prehensive description of the pipeline is given by ^NakosJ ( 2007 1 
and in this paper we provide a brief description of the data pro- 
cessing sequence. 

1 . The sky background, residual features and electronic signa- 
tures were removed from individual (dithered) frames em- 
ploying a combination of comparison frames from which ob- 
jects had been masked. 

2. Individual frames were corrected to a uniform photomet- 
ric sensitivity. Corrections are required to (a) bring all four 
WIRC detectors to the same sensitivity level, (b) remove dif- 
ferential airmass effects, and (c) calibrate all frames to a stan- 
dard photometric system. For this purpose, six photomet- 
ric standard stars (Perss on et al. 1998) were observed dur- 
ing the night of October 26. For the two photometric nights 
these observations were employed to calculate a zero point 
and extinction coefficient. For the non-photometric nights, 
the 2MASS photometry (Skrutskie et al. 2006) was used to 
set the zero point that would shift the instrumental /Tj-band 
magnitudes to calibrated /Tj-band photometry. Full details 
about the calibration can be found in Nakos] p^007| ). 

3. An astro me trie solution was determined for the dithered 
frames and applied to generate a mosaic image. This was 
achieved by comparing the coordinates of point-like sources 
present in the K^, 2MASS and z'-band frames. Having deter- 
mined the astrometric solution, individual frames were com- 
bined to generate four mosaic images. 

4. The source completeness and false detection rate in the 
mosiac images were computed using simulations to deter- 
mine a suitable magnitude limit in the Ks-hand catalog. As 
the Kg catalog is later matched to the R and z' photometric 
catalogs it is the /Tj-band 's faintest magnitude that deter- 
mines the Rz'Ks survey's sensitivity limit. 

Analysis of the mosaiced images indicated that the impact 
of the reduction procedure on the photometric uncertainty was 
at the level of a few hundredths of a magnitude. Exposure mo- 
saics indicating the exact exposure time of pixels in the asso- 
ciated science mosaic were also produced. These frames were 
used to convert the instrumental magnitudes measured on the 
science frames into magnitudes in the Vega photometric system. 
An exposure map of a sub-region of the LCO mosaic image is 
presented in Fig.[T] The sensitivity ratio between the deepest ex- 
posure maps (white) and the shallowest (dark-colored) is on the 
order of three. 



2.4. Constructing ttie K^-band catalog 

Figure [T] indicates that the total exposure time, and therefore 
the background noise, is not constant over the extent of a given 
mosaic. Therefore understanding the spatial noise properties of 
the mosaic is essential for generating the /r5-band object cata- 
log. The limiting source brightness is determined considering: 
(a) the probability of detecting a real source, (b) the number of 




' IRAF is distributed by the National Optical Astronomy 
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Fig. 1. Composite exposure map, generated by combining about 
75 individual exposure maps in which the photometric calibra- 
tion effects (zero point, airmass) have also been incorporated. 
By dividing the science mosaic by its corresponding composite 
exposure map we directly derive the calibrated magnitude of the 
sources in the field. The "sensitivity" of the lighter areas is three 
times higher to that of the dark areas. Sources lying in the over- 
lapping regions of the mosaic (formed by adjacent pointings) 
have a better signal-to-noise ratio. 



false detections as a function of source brightness. The first con- 
dition is required to determine the percentage of sources which 
are missed as a function of the limiting magnitude of the sur- 
vey, while the second condition determines the balance between 
a faint brightness threshold, which will introduce many spurious 
detections, and a bright threshold which will exclude a signifi- 
cant number of interesting objects. 

This study was performed using a code, kindly provided to 
us by E. Labbe (PUC) and F. Courbin (Ecole Polytechnique 
Federale de Lausanne, Switzerland), that simulates the presence 
of unresolved sources in our data. Artificial sources are intro- 
duced into each science mosaic in a regular grid of locations 
modulo a small random offset that varies with each simulation. 
The typical sky area associated with each grid location is 10". 
Source extraction is performed on each simulated image and the 
relative frequency with which the artificial sources are recov- 
ered is stored as the detection probability. This process is then 
repeated for successively fainter artificial sources. The results of 
the simulations, run on all four mosaics, were averaged in order 
to obtain a mean value for the detection probability and for the 
false detection rate. Based on these averaged results, presented 
in Fig. |2j we concluded that excluding detections fainter than 
Ks = 18 (within a Kron-type aperture) from the final /Ts-band 
catalog provides a good balance between a low percentage of 
false detections (on the order of a few percent) and a satisfactory 
completeness of the Ks catalog (80%). In addition to a magni- 
tude cut, sources with erroneous photometry were removed from 
the catalog using a combination of SExtractor flag informa- 
tion and visual inspection. Such rejected sources were typically: 
(a) sources suffering from blending effects due to the presence 
of a counterpart within a few arcseconds, (b) sources close to 
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Fig. 2. Left panel: Average detection probability, over the four science mosaics, as a function of the Ks magnitude (see text for more 
details). Right panel: Average false detection rate as a function of the Ks magnitude. The dashed lines indicate the completeness 
(left panel) and expected percentage of false detections (right panel) for a catalog cut at A'j - 18. 



saturated objects, (c) sources close to the mosaic borders, or (d) 
sources affected by bad columns or other artifacts. 

The final /Tj-band catalog contains the position (RA,Dec) 
and the fixed aperture (3" diameter) photometry of some ~4500 
entries. These entries were selected according to: (a) their Kron 
magnitude (the one used during the simulations), and (b) the de- 
tection probability corresponding to the position of each detec- 
tion in the mosaic images (as derived from the simulations). 

2.5. Matching tlie Rz'K^ catalogs 

Source extraction and photometry was performed using 
SExtractor for the optical R and z' image mosaics using sim- 
ilar parameters and checks as used for the data. Source pho- 
tometry in each band was combined into a single catalog by 
matching the source astrometry in each band within a specified 
2" tolerance. The tolerance value was specified considering the 
astrometric precision of the catalogs and the object density in 
the LCO field: R, z' counterparts were identified for all /Tj-band 
detections, even those with a poorer astrometric accuracy than 
average (0.5"), with a minimum risk of misidentification among 
the few cases of blended objects that were not removed through- 
out the "cleaning" process described above. The final matched 
catalog contains some 3400 sources and general catalog proper- 
ties are presented in Table [T] 

Matched sources were classified as point-like or extended on 
the basis of the R- and z'-band magnitude and the SExtractor 
stellarity index. It should be noted that although the stellarity 
index provides reliable information concerning the morphology 
of relatively bright sources (typical magnitudes R, z' < 22), it is 
less relaible when considering fainter sources. A combination 
of the stellarity index and source magnitude provides a more 
reliable morphological classification. Additional details can be 
found in |Andreon et al.|p004| l. Out of the 3400 sources in the 
matched Rz'Ks catalog some 2450 (~ 72%) were classified as 
extended, while the remainder were classified as point-like. 

The final matched catlogue features 3" aperture photome- 
try for all bands where the photometry is computed using the 
magnitude of Vega as a zero-point. Fig. |3] shows the {R - z') 
vs (z' - Ks) color-color diagram, where the majority of point- 
like and extended sources occupy distinct regions. Despite the 
clear separation of the two populations a small fraction of point 



Table 1. General properties of the /T^-band selected Rz'Ks sur- 
vey. The astrometric and photometric precision are quoted at 1 cr. 



Filter 


Astr. prec. 


Phot. prec. 


Mag limit 


R 


0.3" 


0.02 


22.5 


z' 


0.3" 


0.02 


21.0 


Ks 


0.5" 


0.07 


18.0 



sources overlap with the extended ones and vice-versa, possibly 
reflecting uncertainties in the photometry or the morphological 
classification as well as intrinsically distinct populations. 

3. Multi-wavelength data within the XMDS field 

This paper is primarily concerned with the properties of a quasar 
sample selected using a variant of the KX method, i.e. employ- 
ing the Rz'Ks colors of optically point-like sources in the XMDS 
field. In order to discuss the physical nature of the /TX-selected 
quasar sample we will use photometric information at other 
wavelengths in order to construct low resolution spectral energy 
distributions (SEDs) for the candidate quasars. In addition, in or- 
der to discuss any possible bias in a /TX-selected quasar sample, 
we will also compare the /TX-selected sample to quasar samples 
selected in the mid-infrared waveband. 

We begin by providing a description of the X-ray data 
products provided by the XMDS itself. In addition, the CFHT 
Legacy Survey (CFHTLS) and the Spitzer Wide-area InfraRed 
Extragalacic survey (SWIRE) are two "legacy"-type surveys that 
overlap the XMM Medium Deep Survey (XMDS) field. We pro- 
vide below a brief description of the data products associated 
with these surveys. Detailed information can be found in |Tajer 
|et al.| ( |2007l l and references therein. 

The XMDS was conducted between July 2001 and January 
2003 using the European Photon Imaging Camera (EPIC) on 
board the XMM-Newton satelHte (Turner et al. 200T] IStriider] 
et al. 2001). It consists of 19 overlapping pointings (of typical 
exposure time 20 - 25 ksec) covering a contiguous area of ~ 2 
deg^. The energy range of the EPIC instruments is between 0. 1 to 
12 keV. A detailed description of the pipeline used for reducing 
the X-ray observations and the catalog properties can be found 
in Chiappetti et al. (2005 i. The X-ray source sample used in the 
present paper was generated from the |Chiappetti et al.| ( |2005| l 
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Fig. 3. {R - z') vs (z' - Ks) diagram for the LCO point-like (as- 
terisks) and extended sources (triangles). The computed colors 
are based on the 3" photometry (Vega magnitudes). The point- 
like sources (mainly stars) form a locus well separated from the 
region occupied by the extended sources (galaxies). For com- 
parison we kept the axis convention of Fig. 2 from Warr en et al.| 
( |2000| ), i.e. the vertical axis increases from top to bottom. 



source list by selecting unresolved sources displaying more than 
20 counts in the B (0.5-2) keV energy band. This generated 370 
X-ray sources in the Rz'Ks catalog area. 

The XMM-LSS survey (of which the XMDS is a sub- 
set) is coincident with the Wide Synoptic component of the 
Canada France Hawaii Telescope Legacy Survey 1 (CFHTLS- 
Wl) (Gwyn| |2007). CFHTLS-Wl u*,g',r',i',z' photometry is 
available for approximately 90% of the area covered by the 
Rz'Ks observations. The magnitude limit of the survey, at 
50% completeness, for a point source observed in seeing of 
FWHM=0.8" and a signal-to-noise ratio of 5 (1.5" aperture pho- 
tometry), is 26.4, 26.6, 25.9, 25.5 and 24.8 for the u*, g', r', i' and 
z' bands respectively (AB magnitu des 

SWIRE ( (Lonsdale et"aL][2003| [2004) 1 is the largest Legacy 
Program performed with the Spitzer Space Observatory (fSurace 
|2004 ). Once again, 90% of the sky area associated with the Rz'Ks 
catalog is covered by SWIRE data. The 5cr sensitivity for the 
3.6, 4.5, 5.8, 8.0 and 24 fi m energy bands are 7 .3, 9.7, 27.5, 32.5 
and 450 juJy, respectively (Lonsdale et al. 2003| l. Throughout this 
work we shall use the IRAC photometry measured within a 3" 
aperture; for the MIPS 24 fim band the fluxes presented are those 
measured using a 7.5" diameter. 



the authors. In the present work we implement a variant of the 
original KX technique, using the {R-z') versus (z' -Ks) diagram 
instead. 



4. 1. Seperating quasars from the stellar locus 

Unresolved sources in the Rz'Ks plane essentially consist of 
stars, quasars and unresolved galaxies. The properties of the stel- 
lar locus on the Rz'Ks plane were computed using the Pickles 
(T998) library of 130 stellar spectra (Fig.|4]i. The locus of stellar 
photometry was then extended along the z' - Ks axis to dis- 
play a half-width of 3cr, where cr is the photometric error on the 
(z' - Ks) color index that dominates the photometric uncertainty 
in the stellar locus (see also Table [T]i. Consideration of the rel- 
ative properties of the SEDs of stars and quasars indicates that 
quasars should be located toward z' - Ks colors redder than the 
stellar locus in Fig.|4] We therefore identify candidate quasars as 
unresolved sources displaying 



(R - z') = 0.6 (z - Ks) - 0.243 for {z - Ks) < 2.0, 



(1) 



{R - z) = 2.877 (z - Ks) - 4.797 for 2.0 < (z' - Ks) < 2.4, (2) 
(R - z') = 0.109 (z' - Ks) + 1.983 for (z' - Ks) > 2.4. (3) 
These selection criteria identify 93 candidate quasars. 
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Fig. 4. Simulated stellar colors (shaded region) and distribution 
of the unresolved sources on the Rz'Ks plane. The dashed line 
indicates the adopted color threshold seperating stellar and non- 
stellar sources given by Equations 1-3. 



4. A /^X-selected quasar sample 

The original KX method by [Warren et al.| ( |2000| l made use of 
the V, J and K bands, noting that quasars with sl V - J color 
similar to that of stars would be redder in J-K, a method equally 
well suited for the detection of reddened quasars according to 

^ http://cfht.hawaii.edu/Science/CFHLS/cfhtlsgoals.html 



4.2. SED analysis of quasar candidates 

In the absence of low- to medium-resolution spectroscopy we 
investigate the nature of the candidate quasar sample via anal- 
ysis of their UV to MIR SEDs. The Rz'Ks catalog was there- 
fore matched to the CFHTLS and SWIRE catalogs in order 
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to reconstruct the SED of the objects. The CFHTLS catalog 
was matched to the Rz'Ks detections using a matching radius 
of 2". Thanks to the large overlap between the two sets of 
observations, CFHTLS u* ,g' ,r' ,i' ,z' counterparts were found 
for 86/93 quasar candidates. Furthermore, 90/93 quasar can- 
didates are matched to SWIRE sources using a tolerance of 
1.5" and show emission in the first two IRAC bands (IRAC 1 
= 3.6/im, IRAC 2 = 4.5 ^lm ( [Fazio et al. 2004 )). At longer 
wavelengths, however, the number decreases significantly with 
~ 30% of the candidates showing emission in the IRAC 3 
(5.8 jum) and IRAC 4 (8.0/iOT) band and only 16% with emission 
in the MIPS lA^im channel ( [Rieke e t al. 2004). Therefore, the 
SED of the quasar candidates will include up to 13 photomet- 
ric points (u',g', r' ,R, i',z',z'„i„, K„ IRAC 1, 2, 3 4 and MIPS 24 
bands), spanning the range 0.4 - 24 lum. 

The investigation of the candidate quasars sample was per- 
formed initially via a visual inspection of their SEDs. Quasar 
SEDs exhibit two unmistakable signatures in the UV/optical/IR 
domain: i) the UV/optical-bump customarily attributed to the ac- 
cretion process; and ii) the minimum (in AF^) or break with sub- 
sequent change of slope (in or Fy), occurring at around 1 
micron (see Hatziminaoglou et al. 2005 ; Rich ards et al.| 2006| , 
where the drop of the accretion bump meets with the T ~ 1500 K 
black-body rise (towards the IR), coiTesponding to the graphite 
grains believed to be one of the two main components of the 
dusty tori surrounding the active nucleus. This visual inspec- 
tion yielded 25 quasars among the 93 candidates. The remaining 
sources were characterised as late type stars and (low redshift) 
galaxies, characterised either as early-type or starburst-type. The 
SED of elliptical galaxies peaks around 2 to 3 micron and then 
decreases monotonically towards redder wavelengths, while that 
of starbursts shows an excess in the IRAC 4 band due to the 
strong PAH emission centered at 7.7 micron ( [Weedman et al. 
2006| l. For the quasars, 22/25 (~ 82%) are "classical" (i.e. blue) 
type-1, and one (~ 6%) is of type-2 (as its flux experiences a sig- 
nificant attenuation at bluer wavelengths). Although, at present, 
there is no clear definition of a red quasar ( Canahzo et al.|2006 



Table 2. Results of the SED inspection for the 93 quasar candi- 
dates. 



and references therein), there is no doubt that there exists a frac- 
tion of quasars where the optical spectrum deviates significantly 
from that of a type-1 object (e.g. [Pierre et al.|2001[|Leipski et aUj 
|2007) . Two entries in our quasar sample belong to this category, 
showing a signature of moderate reddening with a flat SED blue- 
ward of the ~1 fj.m inflection. The lack of spectroscopic infor- 
mation and the small number of objects do not allow any firm 
conclusion regarding the origin of the flattening of their optical 
SED. A possible explanation could be the presence of a moder- 
ate amount of dust along the line of sight to the active nucleus. 
In this scenario the dust would affect differentially the blue and 
red light emitted by the quasar, making the SED appear flatter 
compared to that of a blue type-1 quasar, but not as steep as that 
of a type-2. However, larger numbers of sources displaying this 
feature would be needed before claiming any robust conclusion. 
The results of the SED inspection are presented in Table |2] and 
examples of visually classified SEDs are shown in Fig.|5] 



SED type 


Number 


QSO Type-1 


22 


QSO Type-1 (reddened) 


2 


QSO Type-2 


1 


Galaxy 


50 


Star 


18 



CO 

in 




Fig. 6. SWIRE color-color diagram for unresolved sources in the 
Rz'Ks catalog. Sources classified as stars are displayed as aster- 
isks and sources classified as quasars (on the basis of Rz'Ks col- 
ors and SED inspection) are plotted as open squares. According 
to Hatziminaoglou et al.'(2005 ) quasars are expected to be found 



in the region deUneated with the long-dashed line. Visual inspec- 
tion of the SED of the three "stars" within the quasar locus con- 
firms the quasar hypothesis. 



the stellar colors are very close to zero in the Vega photomet- 
ric system. Following their approach, we constructed the MIR 
color-color diagram for all unresolved sources in the Rz'Ks cata- 
log (Fig.|6]). Three sources classified as stars on the basis of their 
Rz'Ks photometry are located in the so-called "quasar locus" and 
visual inspection of their SEDs confirms the quasar hypothesis. 
Of these sources, two are located very close (about 0.2 mag in 
each direction) to the stellar locus on the Rz'Ks plane. The third 
source is an isolated quasar located in the stellar locus of Fig. [3] 
for which we possess a spectroscopic identification (see sec- 
tion 



4.4 1. We finally note that two sources classified as quasars 
on the basis of their Rz'Ks photometry and SED analysis lie out- 
with yet close to the quasar locus defined by |Hatziminaoglou| 
|etaL] ( |2005] l. 



4.3. The MIR color-color diagram 

The KX method seperates quasars and stars on the basis of their 
optical and near-infrared colors. An alternative method applies 
the same test on the mid-infrared colour plane. jHatziminaoglou 
|et al.| (|2005) demonstrated that it is possible to isolate type-1 
quasars using the first three IRAC bands (3.6, 4.5 and 5.8yL(m). 
These bands seperate quasars from stars as they sample the stel- 
lar Rayleigh- Jeans tail of the blackbody spectrum. As a result. 



4.4. Spectroscopy of quasar candidates 

To date it has not been possible to perform a systematic spectro- 
scopic investigation of the KX-selected quasar sample presented 
in this paper. However, a number of spectroscopic programs have 
been performed which overlap the XMDS area with the result 
that spectroscopic information is available for a subset of sources 
presented in this paper. As part of the spectroscopic identifi- 
cation program for the XMM-Newton Survey Science Centre 
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of 1.513, 0.710 and 1.679, respectively. Lower Panel: SED of a star (left), an elliptical galaxy at redshift z - 0.545 (middle) and a 
galaxy showing starburst activity at z =0.127. 



medium sensitivity serendipitous survey Tedds et al. ( 2006 1, the 
two degree field multi object spectrograph (2dF, Lewis et al. 
2002] l targeted some 800 X-ray detected sources with r^s < 22 
including the XMM-LSS area. More details about the 2dF pro- 
gram will be given in a forthcoming X-ray Wide Area Survey 
(XWAS) catalog paper (Tedds et al. in prep.). A total of 22 
sources with 2dF spectra were associated with sources in the 
Rz'Ks catalog. These included one star, 9 quasars and 12 galaxies 
as classified by their Rz'Ks photometry. In each case the spectro- 
scopic classification matched the Rz'Ks classification. A further 
five spectra were available for sources with Ks > 18 and there- 
fore not included in the Rz'Ks catalog. These sources were clas- 
sified as type-l quasars. They are not used further in the analysis, 
except as part of a spectral training set for the photometric red- 
shift analysis (Section |43] l. Three examples of good-quality 2dF 
spectra are presented in Fig.|7] 

A further four sources classified as quasars on the basis of 
their Rz'Ks photometry and SED analysis have been confirmed 
spectroscopically as type-l quasars (source numbers 2, 3, 20 
and 27 in Table 2). Observations were conducted as part of the 
VLT/VIMOS program (080.A-0852B) to investigate X-ray se- 
lected AGN in the XMM-LSS survey ai-ea (Garcet et al. in prep.). 

Therefore, for the ^TX-selected quasar sample, 12/25 sources 
classified as quasars on the basis of Rz'Ks photometry and SED 
analysis have been confirmed spectroscopically. However, it is 
important to note that the spectroscopic programs that observed 



these sources were constructed to investigate X-ray selected 
source samples. It is therefore important at this stage to consider 
the X-ray properties of the /TX-selected quasar sample. 



4.5. X-ray properties of the quasar sample 

As noted in Section|3] there are 370 X-ray sources in the overlap- 
ping XMDS and Rz'Ks catalog area. Of these, 81 X-ray sources 
were matched to a single Rz'Ks source within a 7" matching ra- 
dius. This radius was chosen on the basis of: (a) the density of 
X-ray detections; (b) the accuracy of the X-ray astrometric solu- 
tion (~ 1"); and (c) XMM-Newton's point spread function (6")- 

The relatively low fraction of matched sources, 81/370, is 
largely due to the shallow depth of the Ks band observations 
compared to the optical observations. Approximately 20% of the 
/TX-selected quasar sample, 21/93, are matched to an X-ray de- 
tection. Of these matched sources, 12/13 are spectroscopically 
confirmed quasars and 21/25 are SED classified quasars (note 
that the SED classified quasars include all of the spectroscopi- 
cally confirmed quasars as well). None of the sources classified 
as stars or galaxies as a result of the SED inspection have an 
X-ray counterpart. 

We intend to discuss the X-ray properties of the KX- 
selected quasar sample in greater detail in a subsequent pa- 
per. However, at this point we note that the X-ray-to-optical 
flux ratios of the SED-classified quasars occupy the interval 
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Fig. 7. Examples of good-quality 2dF spectra used to investigate 
Rz'Ks sources. The identification of each object is presented in 
the upper left part of each spectrum; some of the spectral features 
employed to classify each source are also shown. Object # 528 is 
a star, object # 735 a quasar at z = 0.76 and object # 243 a galaxy 
at z = 0.263 ( IGarcet et al.|2007l ). 



0.1 < logiFx/Fu ) < 10.0, typically occupied by AGN ( [Delia 
Ceca et aT]|2004 i. Figure |8] indicates that the four sources clas- 
sified as quasars yet lacking an X-ray counterpart are consistent 
with being drawn from the same population as the X-ray de- 
tected quasars if we consider that application of the above X- 
ray-to-optical flux ratio would place the X-ray fluxes for these 
four sources below the nominal flux limit of the XMDS X-ray 
source catalog. 

Therefore, for the sample of 93 /TX-selected quasars, 21/25 
of the sources classified as quasars from the SED analysis show 
X-ray emission. The remaining four sources are consistent with 
being drawn from the same population. None of the 68 remain- 
ing sources classified as either stars or galaxies from the SED 
analysis is detected in X-rays. We therefore conclude that the 
X-ray properties of the .^TX-selected quasar sample support the 
previous conclusions on the nature of individual sources derived 
from the SED analysis. 

4.6. Photometric redshift analysis 

In the absence of complete spectroscopic informtion for the 
93 sources in the KX-selected quasar sample we employed the 
HyperZ code ( [Bolzonella et al.||2000| l to compute photometric 
redshifts. The spectral templates applied to the data are described 
1 Polletta et al. ( 2007| l and cover the wavelength interval 1 00 nm 



- 100 //m. We tested the performance of HyperZ by first run- 
ning the program on the 14 quasars with 2dF spectra, using vari- 
ous combinations of filters, reddening laws and templates. Based 
on the limiting magnitude of the Rz'K, survey and the maxi- 
mum spectroscopic redshift in our catalog, we set Zma.x - 3.0. 
As HyperZ is limited to 15 templates per run, we selected five 
quasar templates (three type-1 and two type-2) while the remain- 
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Fig. 8. Plot of the X-ray to optical flux ratio for SED and spec- 
troscopically classified quasars and spectroscopically classified 
galaxies (from the 2dF sample) versus X-ray hardness ratio. 
Hardness ratio is computed according to Equation 2 of Garcet 
et al. (2007). The four quasar candidates that are undetected in 
X-rays are plotted as upper limits on F2-10/FR with the median 
HR value of the X-ray detected quasar sample. The region en- 
closed by the dashed rectangle indicates the anticipated locus of 
X-ray emitting stars. 



der were selected among the templates for elliptical, starburst, 
Seyfert and spiral galaxies. 

The performance of the photometric redshift method was es- 
timated on the basis of two parameters: the fractional error Az 

( |2007] l. 



and the Icr dispersion cr, as defined by Tajer et al. 



Zphot 



1 +z 



spec 



] and er'^^-y h""' 



(4) 



where is the number of sources with spectroscopic redshifts. 
After various trials we obtained a photometric redshift for 12 of 
the 14 sources with |Az| < 0.3 and cr, = 0.19. 

The remaining 16 quasars without spectra consists of 15 
type-1 and one type-2 object (#33), characterized as such on the 
basis of their SED. The best-fit template matched by HyperZ 
was in very good agreement with the SED-based classification, 
fitting the ID #33 quasar with a type-2 spectral template and 13 
out of 15 type-1 objects with a type-1 spectral template. Only 
two type-1 objects were fit by a Seyfert 1.8 template (#27 and 
#57). Figure|9]shows the template fitting results for the full sam- 
ple of 93 quasar candidates. Although the upper limit for the 
Ay extinction coeflicient was set to 0.5 no object was found to 
have Ay > 0.4, with 13 out of 19 objects having Ay < 0.2. 
Finally, the B-band absolute magnitude was set to vary between 
-23.0 < Mb < -28.8, but no object was found to have an abso- 
lute magnitude brighter than Mg - -26.22. The final photomet- 
ric redshifts for the 16 newly discovered quasars are presented 
in Table |2] 
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Fig. 9. Left panel: Distribution in the color-color diagram of the initial sample of 93 quasar candidates. The color-coding is based 
on the best-fitting SED returned via the photometric redshift analysis. The dashed line indicates the border between the stellar and 
non-stellar colors used in our selection procedure. Right panel: Color-color diagram containing the quasars detected using jointly 
several methods (2dF, KX, SWIRE). The A'X-selected quasars are marked with filled squares while the open circles denote the 6 
quasars with 2dF spectra not selected by the KX. Three type-1 quasars, selected on the basis of their IRAC colors, are marked as 
open squares. The isolated object at {R - z') ~ 1.7 is the type-2 /TX-selected quasar found in our survey. 



5. Conclusions 

The distribution in Rz'K^ for the 93 quasar candidates is shown 
in the left panel of Fig. |9] The distribution of all quasars found 
in the Rz'Kg catalog, using /TX-selection plus SED analysis, 2dF 
spectra (subsection |4.4[ ) and IRAC colors (subsection |4.3| l is pre- 
sented in the right panel of the same figure. The first conclusion 
to be drawn from these two plots is that the majority of the quasar 
population (type-1) is concentrated in the bluer part of the (R-z') 
axis, with the {R - z') color ranging from 0.2 to 0.8 magnitudes. 
This occurs as the quasar power-law SED produces an excess in 
the /Tj-band, with the result that quasars appear redder than stars 
along the z' - K^ axis. In addition, the same power-law SED 
causes quasars to appear bluer than galaxies. We note the pres- 
ence of a single type-2 quasar detected by the KX-vatthoA. This 
object is much redder, with R - z'~ 1.7, and would have been 
missed if stricter color criteria had been applied, e.g. a bluer cut 
mR-z'. 

To assess the performance of the /TX-method in a 
quantitative manner we adopt the formalism presented 
in IHatziminaoglou et al.| ( [2000) . Defining Nc as the number of 
quasar selected candidates, Nf as the number of quasars con- 
firmed from the candidates and Ne as the number of expected 
quasars (based upon model predictions), one may define the ef- 
ficiency, E, of the technique as 



W2 



(5) 



The KX-msihoA identified 25 quasars in an area of 0.68 deg2 at 
K, < 18. The Rz'K, catalog is 80% complete at K, = 18 and 
therefore the effective surface de nsity of A"X-s elected quasars at 
K, < 18 is 46 + 9deg ^ M addox & Hewett| (2006 ) use model 
predictions to estimate a quasar surface density at Ks < 18 of 
50deg"^. Employing these values, the completeness of the pre- 
selected sample is on the order of 90%. The identification of 25 
quasars from the 93 candidates indicates a confirmation rate on 
the order of 27%. Therefore, employing Equation|5] the overall 
efficiency of the /TX-techniqe as applied in this paper is approx- 
imately 25%. 



Detailed studies of quasar populations represent a key com- 
ponent of the current generation of wide-field optical to MIR 
imaging surveys. Selection by color provides an important 
technique for compiling large samples of quasar candidates. 
However, the application of color criteria alone are unlikely to 
generate uncontaminated quasar samples. SED classification of 
quasar candidates may provide an important technique for pri- 
oritizing the spectroscopic observations of large quasar samples 
- with the potential to improve considerably the efficiency of 
follow-up programs. Throughout the present analysis we have 
applied a visual classification of candidate quasar SEDs. Clearly 
an automated approach will be required for larger samples. 
Large surface IR surveys, such as UKIDSS (Hewett et al. 2006[ l, 
will base their quasar selection on the KX technique. Because 
of the depth of the surveys and the large area coverage (several 
thousand square degrees) it would be advantageous to minimize 
contamination prior to spectroscopic confirmation. 
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Appendix A: Details concerning tlie photo-z 
calculations 

In this appendix we present in detail the parameters used for 
computing the photometric redshifts of the quasars found in our 
survey, for which no spectroscopic information was available. 

HyperZ was fed with 15 templates, taken from the |Polletta 
et al. ( |2007| l library. They were chosen from seven different cat- 
egories of active (AGN) or non-active (ellipticals, spirals) galax- 
ies. Concerning the starburst galaxies, the templates are those of 
Arp 220, M 82 and NGC 6090. The exact number of templates 



per category is presented in Table A. 1 

The minimum and maximum redshifts were 0.7 and 3.0, re- 
spectively, with a step of Az = 0.01. 

In the photo-z computation, the possibility of dust reddening 



was incorporated via the Calzetti et al. (2000 1 reddening law. 
The minimum and maximum reddening values were 0.0 and 0.5, 
respectively, with a step of 0. 1 . 

Constraints were also set on the absolute M^-band magni- 
tude of the sources. We used the gr'-band filter as a reference and 
computed the approximate limits in this band considering the 



formula by [Smith etaLl ( |2002| ): 



B = 3' -H 0.47 {g -r') + Q.\l 



(A.l) 



with -28.0 < Mb < -22.0 and using an average {g' - r')AB 
color index of 0.65. Based on the output of HyperZ, the bright- 
est absolute magnitude was found to be Mg = -25.42 and the 
faintest Mb - -22.5, respectively. 

The cosmological parameters used in this analysis are Hq = 
70kms-'Mpc-',QA = 0.7, = 0.3. 



Table A.l. Distribution of the template categories used for the 
photometric redshift calculations. 



Template type 


Nr 


Type-1 quasar 


3 


Type-2 quasar 


2 


Seyfert 


2 


Starbursts 


3 


EUipt. 


3 


Spirals 


2 
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